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Abstract 


Euler's theorem states that any sequence of finite rotations 
of a rigid body can be described as a single rotation of the body 
about a fixed axis in three dimensional Euclidean space. Th 
usual statement of the theorem in the literature cannot 
extended to Euclidean spaces of other dimensions. Equivalent 
formulations of the theorem are given in this paper 
a way which does not limit them to the three dimensional 
Euclidean space. Thus, the equivalent theorems hold mother 
dimensions. The proof of one formulation presents an atri ^ that 
which shows how to compute an angular-difference ma 
represents a single rotation which is equivalent to the sequenc 
of rotations that have generated the final n-D orientation. T 
algorithm results also in a constant angular-velocity which , when 
applied to the initial orientation, yields eventually the fl 
orientation regardless of what angular velocity generated the 
latter. Finally, the extension of the theorem is demonstrated in 
a four dimensional numerical example. 
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I. EUT.FR '.q THEOREM 


In 1775 
bodies [l] 


Euler published a paper on the rotation of rigid 
In that paper, which was written in Latin, he 


Theor ema. 

Quomodocunquc fphaera circa centrum fuum comer - 
tatur , femper ajjignari pot eft diameter , cuius* 
diredtio in fitu. translato conueniat cum fitiL 
initialL 


Fig. 1 : 


Euler's theorem on the rotation of a rigid 
body as it appeared in the 1775 publication. 


presented the theorem whose photograph is shown in fig. i. The 
theorem states the following; 


In whatever way a sphere is rotated about its center it is 
brings Ppssibie to reckon a diameter about which a rotation 
location^ Sphere lnto coincidence with its original 


A modern formulation of this theorem states [2]: 


A body set of axes at any time t can always 
a single rotation of the initial set of axes 


be obtained 


by 


We prefer to formulate this theorem as follows: 


Regardless of the way a coordinate system is rotated from 

£"r ginal °. rientatlon ' it is always possible to find a 

f* J®*? i X1S ^ - ln space about which a single rotation of the 
initial coordinates ends at the final orientation. 
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Euler' s theorem serves as a cornerstone in attitude 
determination [3-9] and tracking [7,9,10]. In particular, if n is 
a unit vector along the axis of rotation and 0 is the angle by 
which the initial coordinate system has to be rotated in order to 
coincide with the final one, then D(n,0), the transformation 
matrix from the initial to the final coordinate system, is given 
by [7-10] 

D(n, 0) = Icos0 + (l-cos0) nn T - sin0[nx] (1) 

where I denotes the identity matrix, T denotes the matrix 
transpose and [nx] denotes the cross product matrix of n. The 
relationship formulated in ( 1 ) can also be expressed as follows 

D(n, 0) = (2) 

where 0 = n0. The rate of change of the vector quantity 0 as a 
function of w, the angular velocity at which the coordinate 
system rotates, is given by [9,10] 


0 (t) = w(t) + “ 0 (t) xw(t) + 


2 - 0 (t) cot [ 0 (t) / 2 ] 
20 (t ) 2 


0 (t)x[ 0 (t)xw(t) ] 
(3) 


II. ALTERNATE FORMULATIONS OF EULER'S THEOREM 

Ancrular-matrix [ discrete) formulation 

Denote the skew-symmetric matrix [0x] by 0; that is, 

6 = [ 0 x] ( 4 ) 
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where the explicit expression for 9 is: 



Equation (2) can be written as 

D(n,0) - e~ e (6) 

Let the initial orientation of a certain coordinate system with 
respect to some reference system be expressed by the attitude 
matrix D Q . Suppose now that this coordinate system is rotated 
from its initial orientation by the sequence of rotations e 2 , 

* * * • fijc* Denote the cross product matrices which correspond to 
^1' &2 • *•*•' Ik ®l' ••••» fik respectively. Then, in view 

°f (6), the attitude matrix that transforms the reference 
coordinate system to the final one, and which expresses the 
orientation of that system, is given by 

Df “ e ^ . e 63 e 01 D 0 ( 7 ) 

On the other hand 

Df « e ^ D 0 ( 8 ) 

However 

®f J* 6^+ .... +©2+0^ (9) 

We realize that the equivalence between (7) and (8) is another 
expression of Euler's theorem. We can, then, formulate Euler's 
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theorem also as follows: 


Regardless of the way a coordinate system is rotated from 
its original orientation, it is always possible to express 
the final orientation of the system by the attitude matrix 
Df where 

-e f 

D f = e D 0 

and ©£ is a skew-symmetric matrix. 

We call this formulation discrete because the rotation expressed 
by © f is equivalent to the discrete k rotations expressed by the 
individual (i=l,2, k) matrices. 

Angular-rate (continuous) formulation 

Euler's theorem gave rise to (1) and (3) which indicate how to 
find the orientation of a coordinate system, at any given time 
t^, with respect to its initial orientation at time t Q if w(t) , 
the history of its rate of rotation, is known for t Q <t<t^; 
namely, w(t) is used in (3) to solve for ,0(t) and then, that 
solution is used in (1) to obtain the required orientation 
specified by D(t f ) . The attitude matrix D(t f ) can be computed 
in yet another way, since the rate of change of D(t) as a 
function jf w(t) is given by the well known matrix differential 
equation 

D(t) = - [w(t)x]D(t) (10) 

The matrix [w(t)x] is defined on the components of w(t) when the 
latter is resolved in the changing (final) coordinate system. We 
also denote this matrix by W(t) ; that is, 
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W(t) = [w(t)x] 


( 11 ) 


Hence (10) can be written also as 

D(t) = - W(t)D(t) (12) 

The explicit expression for W(t) ( or [w(t)x] ) is given by the 
skew-symmetric matrix: 

0 -w 3 (t) w 2 (t)~ 

W(t) = w 3 (t) o -w^t) ( 13 ) 

-w 2 (t) w x (t) 0 

There is, then, an equivalence between the pair (1) and (3) on 
the one hand, and (12) on the other hand. 

Euler's theorem states, basically, that there always exists a 
vector 0(t f ) which specifies the orientation regardless of 
which w(t) generated that ©(t f ) . Consequently, any w(t) which 
satisfies the following two conditions, rotates the initial 
coordinate system into the same orientation 

w(t) _ * 

w ( t ) n (14. a) 

f tf 

j w(t) dt = 0 (tf) (14. b) 

^•o 

While the truth of the last proposition is self evident, it can 
also be easily verified by solving (3) for any angular rate w(t) 
which satisfies conditions (14) . Since any w(t) which satisfies 
(^-4) rotates the initial coordinates into the same orientation. 
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then certainly the constant angular rate vector specified by 


w 


e(t f ) 


(15) 


rotates the initial coordinates into the same orientation and 
since according to Euler's theorem such 0(tf) always exists, then 
such constant w exists too. Finally, since such a constant 
exists then, following (4) and (11), there also exists a 
corresponding constant matrix, W, 


W = ----- © f (16) 

'"f t o 

which when used in solving (12), yields the D(t^) that 
corresponds to 0(tf) . Therefore, in view of the equivalence 
between the pair (1) and (3) on the one hand, and (12) on the 
other hand, we can phrase an equivalent formulation of Euler's 
theorem as follows: 

Regardless of what matrix W(t) generated D(tf), it is 
always possible to find a constant matrix W which generates 
the same D(tf) . 

We call this formulation continuous because it relates to the 
continuous change of the orientation as a result of the existence 
of an angular rate at which the orientation changes. 


III. REPRESENTATION OF ROTATIONS IN n-D 

Denote the dimension of an Euclidean space by n. The rotation 
matrix in n-D, being a square matrix, consists of n 2 elements. 


However, the orthogonality of the matrix imposes (n+l)n/2 
constraints on it. Consequently a rotation matrix in n-D has only 
m=(n-l)n/2 independent parameters. That is, a rotation matrix in 
n-D is defined by exactly m=(n-l)n/2 parameters. Consider now the 
3-D rotation. As indicated by (1) and indeed as stated by Euler 
[1,11] the vector 9(t f ) in its three components contains the 
necessary and sufficient information for specifying the 3-D 
rotation. Similarly the orthogonal rotation matrix D(t f ) contains 
three independent parameters although it has nine elements. So 
the 3-D case is unique in that n=m and the rotation can be 
described by either a vector or a matrix. In all other 
dimensions, though, n^m and since m parameters are needed to 
define the rotation, a vector with its n elements cannot 
specify a rotation. Rotation matrices though, with their m 
independent parameters, do specify the rotation. As a 
consequence of this discussion, it is concluded that the 
original version of Euler's theorem or any of its variants 
presented in Section I are not extendible to n-D while the 
alternate formulation of Euler's theorem given in the preceding 
section may be extended to n-D. 


IV. EULER'S THEOREM IN n-D 

In view of the conclusion drawn in the last section, the 
general formulations of Euler's theorem in n-D are that given in 
Section II. Let us first address the angular-matrix (discrete) 
formulation and rephrase it in a more general frame by the 
following theorem: 
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Theorem 1: Given the arbitrary unitary matrices D Q =D(t Q ) and 

Df=D(t f ), then D f can always be expressed In the 
form: 

D f - e" e D q (17) 

where 0 is a skew Hermit lan-matrix. 

Since the rotation matrix is orthogonal, the angular-matrix 
formulation of Euler's theorem is a special case of this theorem. 
Although the theorem is not new [see e.g. exercise 4 on p. 346 of 
ref. 12], for the sake of completeness, we present here a proof 
of the theorem. 

Proof : Define the unitary matrix 

D = D f Dj (18) 

where + denotes the conjugate transpose of a matrix. Since D is 
unitary, it is also normal and as such it has n orthogonal 
eigenvectors (see theorem A1 in the Appendix) . Define a matrix V 
whose columns are the eigenvectors of D. Then V is unitary. Since 
the eigenvectors of D form an orthonormal set, then 


D = V G V + (19) 

where G is the diagonal matrix of the eigenvalues of D (see A2) . 

Now since D is unitary, its n eigenvalues g 1 ,g 2 , /9 n 

lie on the unit circle of the complex plane (see A3) ; that is, 

g^ ~ e^i 1-1,2, ... , n 


where j-(-i) 1 / 2 and is the phase of the i th eigenvalue. Let us 
form a diagonal matrix S 
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® “ diag{-0 lf -0 2 , ,-0 n } 


( 20 ) 


Then, obviously. 


G = e"3 ffi 


( 21 ) 


Next we define a constant matrix e as 


© = V j® V + 


( 22 ) 


then (see A4) 


r® = e" V 3® V + _ v v+ 


(23) 


Substituting (21) into (23) we obtain 


e“® = V G V + 


(24) 


A comparison between (24) and (19) yields 


D = e~® 


(25) 


Then from (18) and (25) we obtain 

D f " e_9 D o 

To complete the proof we still have to show that 0 is skew- 
Hermitian. From (25) 


D + = e- 9+ 


(26) 


Also (see A5) 


D -1 = e® 


(27) 


but, since D is unitary 
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therefore the right hand side of (26) is equal to that of (27) , 
consequently 

e + = -© 


This completes the proof. 

With Theorem 1 on-hand we are ready now to address the 
angular— rate (continuous) formulation of Euler's theorem. Here 
too we rephrase the latter in a more general frame as follows: 

Theorem 2: Given the arbitrary unitary matrices D =D(tq) and 

D f =D(t f ) , then D f can always be obtained as a 
solution of (12) with the initial condition D Q 
where W is a constant skew-Hermit ian matrix. 

Note that due to the orthogonality of the rotation matrix, the 
angular-rate formulation of Euler's theorem constitutes a special 
case of the last theorem. The following proof of the latter is 
based on the former theorem. 


Proof : From Theorem 1 


D f = e “ 9 D o 


On the other hand, the solution (12) when W is constant is 


-W(t 0 -tf) 

D f = e v ° D c 


the equality of these yields 


W = 


tf-t G 


(28) 


that is; no matter what W(t) generated Df, we can always find a 
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constant W according to (28) for which the solution of (12) with 
the initial condition D Q yields D f at t f . Since e is skew- 
Hermitian then, from (28), W is skew-Hermitian too. This 
concludes the proof. 

While we can consider W as an average velocity computed by 
(28) , we note however that W is not the average of W(t) ; that is, 

w t Z--Z- f W(t) dt (29) 
t f" t o t J 0 

This inequality is known in 3-D as non-commutativity. (Another 
expression for the non— commutativity of rotations is the 
inequality expressed in (9) ) . 

The above theorems extend Euler's theorem in two ways. First 
they deal with the general n-D rather than the 3-D case and 
secondly they extend Euler's theorem to the unitary (complex) 
transformation. Euler's original formulation is, then, a special 
case of the above theorems. 

V. NUMERICAL EXAMPLE 

To demonstrate the facts pointed out in the preceding section 
we bring a fourth order example in which we show how D f which is 
®^tained as a result of the solution of (5) for a certain time 
varying angular velocity matrix, W(t) , can be obtained by the 
solution of (5) with a constant angular velocity matrix, W. 

For simplicity we deal with a special unitary matrix; namely, 
with an orthogonal one. Also for simplicity and with no loss of 
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generality, we choose D 0 =I. We first use the following time 
varying W(t) to solve (12) from t o =0. to tf=0.5 sec 


W(t) = 


" 0 . 

1 . 5t 1 . 5t 2 

0 . 8t 

-1 . 5t 

0. -0. 9sin (6 . 28t) 

-0. 95/(1 

-1 . 5t 2 

0.9sin(6.28t) 0. 

0.75 

-0 . 8t 3 

0. 95/(1. -t) -0.75 

0 . 


3 “I 


(30) 


this yields the following solution at t^: 


0.98130682 

0.18388549 

0.04691911 

-0.03196326 


-0.15805594 

0.76180341 

-0.10221727 

-0.61985926 


-0.08266215 

0.21777062 

0.96379421 

0.12978307 


-0.07226489 

0.58173674 

-0.24176631 

0.77324588 


The following eigenvector matrix of Df is obtained using the 
EIGZF routine of the IMSL library: 


V 


-0 . 038+ j 0 . 161 
0 . 683+ j 0 . 

-0 . 17 0+ j 0.184 
_0 . 056+ j 0 . 664 


-0 . 038- j 0 . 161 
0 . 683+ j 0 . 

-0. 170-j0. 184 
0 . 056- j 0 . 664 


0 . 687+ j 0 . 
0.037-j0. 177 
-0.052-j0.659 
-0 . 151+ j 0.181 


0. 687+jO. 
0.037+j0. 177 
-0 . 052+ j 0 . 659 
-0. 151-jO. 181 


The same routine also yields the following eigenvalues of D f : 


g x = 0 . 7452+j 0 . 
g 2 = 0. 7452-jO. 
g 3 = 0 . 9949+ j 0 . 
g 4 = 0. 9949-jO. 


6668 

= e jO. 7300 

6668 

= e -j0.7300 

1012 

= e jO. 1013 

1012 

= e ~ j 0 . 1013 


The diagonal matrix j® is computed according to (20) and used in 
(22) and (28) to compute the following constant angular velocity 
matrix: 
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w 


0.00000000 

-0.37147707 

-0.12387442 

-0.04981448 


0.37147707 

0.00000000 

0.35388342 

1.31771209 


0.12387442 

-0.35388342 

0.00000000 

-0.40459905 


0.04981448 

-1.31771209 

0.40459905 

0.00000000 


(31) 


When now (12) is solved with the initial condition D Q =I starting 
at to 11 " 0 * then, as expected, D f is the solution obtained at time 
t=0.5. Finally, to demonstrate the inequality presented in (29) 
we average W(t) given in (30). The resulting average is: 


1 /tf 

0 . 

0.375 

0.125 

0.025 

W(t) dt = 

-0.375 

-0.125 

0 . 

0.573 

-0.573 

0 . 

-1.317 

0.750 

** ^ 

-0.025 

1.317 

-0.750 

0 . 


Obviously the latter matrix differs from the constant angular 
matrix, W, given in (31) . 


V. CONCLUSIONS 

Euler's fundamental theorem on the ability to describe any 
orientation of a rigid body as a single rotation and the various 
known versions of this theorem cannot be directly extended to 
other dimensions because all known formulations hinge on the 
concept of axis of rotation which does not exist in dimensions 
other than three. Nevertheless, when it is recognized that the 
general n— D rotation is characterized not by an axis of rotation 
but rather by an angular— difference matrix or by an angular 
velocity matrix, Euler's theorem can be reformulated in 3-D in 
ways which are equivalent to the other known formulations and 
then the new formulations can be extended to n-D. In this work we 
presented the new formulations in 3-D and then we proved that 
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they hold for any dimension. One of the new formulations states 
that no matter what was the sequence of rotations that resulted 
in the final orientation, it is always possible to express the 
rotation matrix as an exponential function of a skew-symmetric 
angular— difference matrix. The other new formulation of Euler's 
theorem states that no matter how the angular velocity matrix 
changes as a function of time, we can always find a constant 
angular-velocity matrix which will result in identical 
orientation change over the same time interval. 

The proof of the theorems supplied the algorithm needed to 
compute the angular-difference matrix and the equivalent constant 
angular-velocity matrix once the initial and the final attitudes 
as well as the time interval are given. To demonstrate the new 
formulations of the theorem and their extendibility to dimensions 
other than 3 we used the algorithm to solve a 4-D example. The 
example clearly demonstrates the ability to reach the same final 
orientation using a constant angular velocity matrix. 


Appendix 

This appendix lists some known theorems which are used in the 
proof of the theorem on the extended Euler theorem. 

Al : A set of n orthonormal eigenvectors can be found for an nxn 

normal matrix. [See p. 76 of Ref. 13]. 

A2 : A matrix can be reduced to a diagonal matrix by a similarity 

transformation if and only if a set of n linearly 
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independent eigenvectors can be found. [See p. 72 of Ref. 

13] . 

A3: The eigenvalues of a unitary matrix have absolute value 1. 

[See p. 129 in Ref. 14]. 

M: If A = TJT -1 , then f (A) = Tf(J)T _1 . [See p. 80 of Ref. 13]. 

A5: If A = e B then A~^ = e -B . 

Proof : Since B and -B commute, then e B e _B = e B-B = I 
hence 

Ae -B =I 

thus 

A -1 =e~ B 
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